Hydrodynamic View of Wave-Packet Interference: Quantum Caves 
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Wave-packet interference is investigated within the complex quantum Hamilton- Jacobi formalism 
using a hydrodynamic description. Quantum interference leads to the formation of the topological 
structure of quantum caves in space-time Argand plots. These caves consist of the vortical and 
stagnation tubes originating from the isosurfaces of the amplitude of the wave function and its 
first derivative. Complex quantum trajectories display counterclockwise helical wrapping around 
the stagnation tubes and hyperbolic deflection near the vortical tubes. The string of alternating 
stagnation and vortical tubes is sufficient to generate divergent trajectories. Moreover, the average 
wrapping time for trajectories and the rotational rate of the nodal line in the complex plane can be 
used to define the lifetime for interference features. 

PACS numbers: 03.65.Nk 



One of the most fundamental but intriguing micro- 
scopic effects is quantum interference, the observable fea- 
ture arising from the coherent superposition of quantum 
probability amplitudes. Quantum interference is involved 
in a very wide range of experiments arising from myriad 
applications. Just to mention some of them, there are 
superconducting quantum interference devices [l[ , coher- 
ent control of chemical reactions Q , atom and molecular 
interferometry [3j (including Bosc-Einstcin condensates 
[3]), or Talbot/Talbot-Lau interferometry with relatively 
heavy particles (e.g., Na atoms [j| and Bose- Einstein con- 
densates @). However, despite all this experimental and 
theoretical work, very little attention beyond the impli- 
cations of the superposition principle has been devoted 
to understanding quantum interference at a more funda- 
mental level 0). 

In this Letter, we focus on the hydrodynamical inter- 
pretation [8, 9] of experiments of this type by introducing 
complex quantum trajectories originating as characteris- 
tics of the solutions of the complex quantum Hamilton- 
Jacobi equation 13, 3, [12, ESQ) Ell- As recently shown 
, interference effects on the real axis may be described 
in terms of the superposition of amplitudes carried by 
approximate (low-order) complex quantum trajectories. 
Unfolding of the dynamics from real space into the com- 
plex plane yields unexpected and surprising features, in- 
cluding what we term quantum caves. These caves are 
topological structures developed around curves in com- 
plex coordinate space where the total wave function and 
its first derivative are zero {nodes and stagnation points, 
respectively). Quantum caves are then displayed in 3D 
Argand plots (the third dimension being time), where 
vortical and stagnation tubes form around nodal and 
stagnation curves (which arise from the time-evolution 
of nodes and stagnation points, respectively), displaying 
analogies to the stalactites and stalagmites of real geolog- 
ical caves. 

The equation of motion for complex quantum trajec- 



tories arises after substituting the complex-valued wave 
function in the form ^(x,t) = exp[iS(x,t)/h] into the 
time-dependent Schrodinger equation. This yields the 
complex- valued quantum Hamilton- Jacobi equation 
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where S(x,t) is the complex action and the last term is 
the complex quantum potential, Q(x,t). For the system 
studied here, no external interaction potential is assumed 
(i.e., V = 0). Quantum trajectories are then developed 
from the guidance condition p(x, t) — dS(x, t) /dx, which 
defines the quantum momentum function (QMF). By an- 
alytical continuation, the x variable is extended to the 
complex plane through the z — x + iy complex vari- 
able (time remains real-valued) and complex quantum 
trajectories are determined from p(z,t) — dS(z,t)/dz = 
(H/i)dln^(z,t)/dz. Two kinds of singularities are espe- 
cially relevant: (i) nodes of the wave function, which cor- 
respond to poles of the QMF, and (ii) stagnation points 
[17| , which occur where the QMF is zero and correspond 
to points where the first derivative of the wave function 
is also zero. In addition, caustics are related to free wave- 
packet propagation [lj|. 

To illustrate the formation of vortical and stagna- 
tion tubes and quantum caves, we consider the head- 
on collision of two one-dimensional Gaussian wave pack- 
ets 0, Us} . Despite its simplicity, this analytical prob- 
lem is a representative of other more complicated, realis- 
tic processes characterized by interference. This process 
can be described by the total wave function, ^(x,t) = 
^L{x,t) + yj R {x,t) (L/R denotes left/right), which is an- 
alytically continued to the complex plane to give ^(z,t). 
Each partial wave is represented by a free Gaussian wave 
packet, 



ip(x,t) = A t i 



-(x—Xt) 2 /4<T t ao+ip(x-x t )/h+iEt/h 



(2) 
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FIG. 1: Quantum caves for head-on collision of two Gaus- 
sian wave packets. These caves are formed with the iso- 
surfaces = 0.053 (pink/lighter gray sheets) and 

\d^{z, t)/dz\ = 0.106 (violet/darker gray sheets). The com- 
plex quantum trajectories launched from two branches of the 
isochrone reach the real axis at t = 5. 

where, for each component, A t = (2naf)~ 1 / 4 and the 
complex time-dependent spreading is given by a t = 
(To(l + iht/2-ma 2 ) with the initial spreading cto- Due to 
the free motion, xt = xq + vt (v = p/m is the propaga- 
tion velocity) and E — p 2 /2m. We also consider the case 
where the relative propagation velocity is larger than the 
wave packet spreading rate, h/2maa 7\. From now on, 
all quantities will be given in atomic units (H= m = 1). 

The following initial conditions are used: x QL = — 10 = 
— Xqr, vl = 2 = — vr and Co = a/2, and maximal 
interference occurs at t = 5 in real space. In Fig. [T] 
complex quantum trajectories together with the isosur- 
faces \$f(z, t)\ — 0.053 (pink/lighter gray sheets) and 
\d^(z,t) / dz\ = 0.106 (violet/darker gray sheets) from 
t = to i = 10 are shown in a 3D Argand plot. 
Around nodes and stagnation points, tubular shapes de- 
velop (pink/lighter gray and violet /darker gray tubes, 
respectively), which alternate with each other and whose 
centers correspond to vortical and stagnation curves, re- 
spectively. The sharp features and well defined vertical 
tubes observed in Fig. [TJ reminiscent of stalactites and 
stalagmites, lead us to call these plots quantum caves. 
Interference leads to the formation of quantum caves and 
produces this topological structure. 

As seen in Fig. [T] and, in more detail, in Fig. [2ja), 
the complex trajectories display counterclockwise helical 
wrapping around the stagnation tubes, while they are hy- 
perbolically deflected or "repelled" when they approach 
the vortical tubes enclosing the QMF poles. This intri- 
cate motion depicts the probability density flow around 
the vortical and stagnation tubes. Trajectories launched 
from different initial positions may wrap around the same 



stagnation curve and remain trapped for a certain time 
interval. As time proceeds, these trajectories separate 
from the stagnation curves in analogy to the decay of a 
resonant state. Therefore, the whole process shows long- 
range correlation among trajectories arising from differ- 
ent starting points. 

The QMF can be viewed as a vector field in the com- 
plex plane, p = p x + ipy, and we can compute its di- 
vergence and vorticity along a complex quantum tra- 
jectory, which describe the local expansion or contrac- 
tion and rotation of the quantum fluid, respectively. By 
the Cauchy-Riemann equations, the first derivative of 
the QMF becomes dp/dz = (T + iQ)/2, where T = 
V • p = dp x /dx + dpy/dy is the divergence of the QMF 
and SI = |V x p \ = (dp y /dx — dp x /dy) is the vorticity. 
Moreover, the complex quantum potential in Eq. |T]) can 
be expressed in terms of divergence and vorticity by 

Q(*,t) = JL^ = JL(r + »n). (3) 

Ami oz Ami 

Figure [2ja) shows trajectories 1, 2, 3 and 4 launched 
from the isochrone which arrive on the real axis at t = 5 
(maximal interference), and Fig. [^b) presents the time 
evolution of the divergence and vorticity of the QMF 
along trajectory 1. When the particle approaches the 
vortical curve at position a, it experiences a repulsive 
force provided by the pole of the QMF and the trajectory 
displays hyperbolic deflection. As shown in Fig. [DJb), T 
and SI display the first sudden spike. Then, this parti- 
cle is trapped by the stagnation curve between two vor- 
tical curves. When the trajectory approaches turning 
points (6, c and d), the particle's velocity undergoes rapid 
changes and this produces sharp fluctuations in T and Q. 
From Eq. ^ , the quantum potential is larger near these 
positions. Finally, as the particle departs from the stag- 
nation curve, it experiences a repulsive force provided by 
the pole and the trajectory displays hyperbolic deflection 
at position e. The whole process indicates important dy- 
namical activity, which is lacking within the real-valued 
version of this problem (where no divergence or vorticity 
can be defined). 

The wrapping time for a specific trajectory can be de- 
fined by the interval between the first and last minimum 
of SI, and the positive vorticity within this time inter- 
val describes the counterclockwise twist of the trajectory. 
The sign of T indicates the local expansion or contrac- 
tion of the quantum fluid when it approaches or leaves a 
turning point, respectively. Within this time interval, the 
particle obviously feels the presence of stagnation points 
and nodes, and the trajectory displays the interference 
dynamics. From Fig. [2f b) , the wrapping process lasts 
from t « 3.7 to t « 6.9. In addition, trajectories 1 and 
2 wrap around the same stagnation curve with differ- 
ent wrapping times and numbers of loops. The wrapping 
time around a stagnation curve is determined by T and SI, 
which are used to characterize the turbulent flow. Thus, 
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FIG. 2: (a) Trajectories 1, 2, 3 and 4 launched from the isochrone which arrive on the real axis at t — 5 display hyperbolic 
deflection around the vortical curves (pink/light gray) and helical wrapping around the stagnation curves (violet/dark gray). 
Trajectories 4 and 5 diverge at position "V" near the vortical curve, (b) Divergence and vorticity of the QMF along trajectory f . 



the average wrapping time for those trajectories reaching 
the real axis at the time of maximal interference can be 
used to define the "life-time" for the interference process 
observed on the real axis. 

Trajectories 2 and 3 start from the isochrone with the 
initial separation Azo ~ 0.3, wrap around different stag- 
nation curves, and then end with the separation at t = 10 
Az » 0.8. These two trajectories avoid the vortical curve 
and this greatly increases the separation between them. 
This behavior is consistent with what one observes when 
looking at the quantum flow in real space: the trajectory 
distribution is sparse near nodes of the wave function 
and dense between two consecutive nodes. In addition, 
trajectories 4 and 5 start with slightly different initial 
positions, Azo = 0.01, and they suddenly separate at 
position "V" near the vortical curve. This leads to the 
continuously increasing separation between them and a 
positive Lyapunov exponent, analogous to the case re- 
ported in real space [13j. The alternating structure for 
the vortical and stagnation tubes (similar to that for the 
nodal point -X-point complex in Bohmian mechanics in 
2D real space) thus leads to divergent trajectories and 
may generate chaos fl9j | . 

Time-dependent nodal positions in the complex plane 
can be determined analytically by solving the equation 
^(z,t) = 0, which renders 



,(*) 



in (n + 1/2) 



[imv/h - (.T - vt) I (2<7oct)] 



(4) 



where n = 0, ±1, ±2, . . . Splitting this expression into its 
real and imaginary parts, z n (t) — x n (t) + iy n (t), yields 
the analytical expression for the angle of the nodal line 
with the positive real axis, 6(t) = tan -1 [y n (t)/x n (£)], 
(which is independent of n), and this describes the time 
evolution of the nodal line. In addition, the time evolu- 
tion of the nth node is given by y n = (2mvo^jhx^)x n — 

(2n+l)(7rog/so). 

Figure [3] shows the time-dependent string of stagna- 
tion points and nodes and nodal trajectories in the com- 



plex plane. At t = 0, these two wave packets are far 
away from each other; however, their tails interfere in 
the complex plane. This contributes to the string of 
stagnation points and nodes, and the initial angle of the 
nodal line (which is perpendicular to nodal trajectories) 
is O = tan~ 1 (—hx /2<mvoQ) = -51.34°. Then, the nodal 
line rotates counterclockwise and crosses the real axis at 
t = 5, where the total wave function displays maximal 
interference. At this time, y n — and we recover the 
expression for the positions of nodes, x n = (n + l/2)A/2 
(A = h/mv). After t = 5, these two wave packets start 
to separate, and the nodal line continues to rotate coun- 
terclockwise away from the real axis. However, these two 
wave packets still interfere with each other in the complex 
plane. When t tends to infinity, the angle of the nodal line 
approaches 9^ = ian^ 1 {2mv(jQ / Kxq) — 38.66° and this 
line becomes parallel to nodal trajectories. The nodal line 
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FIG. 3: Evolution of the nodal line at t = 0, 2.5, 5, 7.5, 10 
(black solid lines) and t — oo (black dashed line); the arrows 
indicate the rotation direction. Nodes and stagnation points 
are denoted by dots; nodal trajectories are shown as dotted 
lines passing through the nodal points. 
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rotates counterclockwise [18| from 9q to a limiting value 
6*00 with an angular displacement AO — doa — 9q = ir/2, 
and its intersections with nodal trajectories determine 
the positions of nodes. In addition, the distance be- 
tween stagnation points and nodes increases with time. 
In particular, if both wave packets are initially very far 
apart (xo — > oo) but move with a finite velocity v, or 
they are separated by an arbitrary finite distance with 
v = 0, then the nodal line will end up aligned with 
the real axis. Interference features are observed on the 
real axis only when the nodal line is near the real axis. 
For the case shown in Fig. [3l this occurs between about 
0(3.52) = -10° and 0(7.32) = 10°, so that the "lifetime" 
for the interference features is about At = 3.8. There- 
fore, the lifetime of interference features observed on the 
real axis is determined by the rotation rate of the nodal 
line in the complex plane, and this rate d9{t)/dt decays 
monotonically to zero as t — > oo. 

The complex quantum trajectory method provides an 
insightful alternative to the traditional analysis of quan- 
tum interference phenomena in real space. In Bohmian 
mechanics, when two or more real coordinates are in- 
volved, quantum vortices form around nodes in the wave 
function and streamlines surrounding the vortex core 



form approximately circular loops [20|, [2l|. The inter- 
ference of two wave packets in one real coordinate leads 
to the formation of nodal structure, but quantum tra- 
jectories close to nodes forming at the maximal inter- 
ference time do not display vortical dynamics 



quantum potential near these nodes forces these trajecto- 
ries to avoid these regions and to exhibit laminar flow in 
space-time plots. In contrast, complex quantum trajecto- 
ries displaying helical wrapping and hyperbolic deflection 
undergo turbulent flow in the complex plane. This coun- 
terclockwise circulation of trajectories launched from dif- 
ferent positions around the same stagnation tubes can be 
viewed as a resonance process in the sense that during in- 
terference some trajectories keep circulating around the 
tubes for finite times and then escape as time progresses. 
On the other hand, in conventional quantum mechan- 
ics, the interference pattern transiently observed on the 
real axis is attributed to constructive and destructive in- 
terference between components of the total wave func- 
tion. In contrast, within the complex quantum trajectory 
formalism, two counter-propagating wave packets are al- 
ways interfering with each other in the complex plane. 
This leads to a persistent pattern of nodes and stagna- 
tion points which is a signature of the "quantum coher- 
ence" demonstrating the connection between both wave 
packets before or after interference fringes are observed 
on the real axis. The interference features observed on 
the real axis are connected to the rotational dynamics of 
the nodal line in the complex plane. Therefore, the aver- 
age wrapping time for trajectories and the rotation rate 
of the nodal line in the complex plane provide two meth- 
ods to define the interference lifetime observed on the 



real axis. This analysis demonstrates that the complex 
quantum trajectory method provides a novel perspective 
and leads to new insights for analyzing and interpreting 
quantum mechanical problems. Finally, similar conclu- 
sions are drawn when the spreading velocity of the wave 
packets is greater than their propagation velocity. 
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